GENERALIZED LIE-BACKLUND THEOREM FOR 
LIE CLASS co = 1 OVERDETERMINED SYSTEMS 

BORIS KRUGLIKOV 

Abstract. In this paper we prove a version of Lie-Backlund the- 
orem for overdetermined systems of scalar PDEs, whose general 
solution depends on 1 function of 1 variable. This generalizes the 
case of involutive system of the second order on the plane treated 
by E.Cartan in 1910. Many examples are provided. 



Introduction 

Consider the space J k n of fc-jets of sections of a vector bundle it : 
E — > M. This manifold carries the canonical Cartan distribution Ck 
(also known as higher contact vector distribution). Diffeomorphisms of 
J k n preserving the Cartan distributions are called Lie transformations. 

The classical Lie-Backlund theorem states that a Lie transformation 
F : J fc 7r —7- J k n is the prolongation of a Lie transformation / of J e 7c: 
F = f( k ~ e \ where e = 1 and / is a contact transformation of J lr K in the 
case rank7r = 1, while e = and / is a diffeomorphism of J°tt (point 
transformation) for rankff > 1. 

Various generalizations of this to Lie transformations of differential 
equations £ C J k n have been discovered since, and they became known 
as Lie-Backlund (type) theorems. Namely an internal transformation is 
a symmetry of the induced distribution = CdTS on S. Lie-Backlund 
theorem holds for £ if the internal symmetries coincide with the exter- 
nal symmetries, which are (restrictions of) those Lie transformations 
of the ambient jet-space that preserve £. 

An instance of such theorem for scalar second order equations was 
proved in IC2I (parabolic equation) and |GK] (hyperbolic and elliptic 
PDE). For many classes of equations (C-general, normal) this phenom- 
enon was established in |KLV| IAKU] . 

But Lie-Backlund theorem does not always hold. It is easily seen to 
fail for a single scalar PDE of 1st order F(x 1 , . . . , x n , u, ui, . . . , u n ) = 
for which the internal symmetry group is infinite-dimensional. An- 
other important counter-example is the Hilbert-Cartan equation y' = 
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(z") 2 , for which the external symmetry group is 9-dimensionaO vs. 14- 
dimensional internal symmetry group, though the latter coincides with 
the space of generalized Lie-Backlund symmetries |AKO| IKer] . 

More general Monge equations = F(x, y, . . . , ?/ m_1 ), z, . . . , z^) 
were considered in [M], If m = l,n > 2 or m > 1 then Lie-Backlund 
theorem holds provided the equation is non-degenerate (this condition 
meaning that F is not affine in z^ 1 cannot be dropped) and the external 
symmetries concern the ambient mixed jet-space J m,n (R, R x R). 

The Monge equations are included into a more general class of Monge 
systems (systems of ODEs with the minimal degree under-determinacy) , 
which are the natural target of reduction for so-called Lie class ou — 1 
compatible overdetermined PDE systems. These are the systems £ 
whose general solution depends on 1 function of 1 variable. 

Lie class u> = 1 PDE systems are integrable by ODE methods [Lj 



K3I , and they often arise as symmetry reductions of more complicated 
PDEs. Lie-Backlund theorem clearly fails for Lie class u> = 1 overde- 
termined PDE systems, since the internal symmetry group is always 
infinite-dimensional while the external group is usually not. 

An important example of such £ constitute involutive overdeter- 
mined systems of 2nd order scalar PDEs on the plane considered by 
E.Cartan ("le theoreme imporant" IC2I , §26]). He indicated a general- 
ization of Lie-Backlund theorem to this case, namely the internal group 
is changed to the symmetry group of the reduction (M, A) of (£,Cs) 
by the Cauchy characteristic n C Cs and then it is bijective with the 
external symmetry group. 

This was generalized to involutive 2nd order PDE systems in n- 



dimensions in W\ Y 2 . The purpose of this paper is to prove 



Theorem . Let £ be a Lie class u = 1 overdetermined compatible 
system of PDEs of orders greater than 1. If £ is sufficiently non- 
degenerate, then the symmetry algebra of the reduction (M, A) is equal 
to the Lie algebra of external symmetries of £ . 

The requirement of sufficient non- degeneracy is technical and will be 
formulated, after we introduce some preliminary material, in Section 
HI Then in Section Owe will prove the main result. 

We will apply it to calculate the group of contact transformations 
of some model overdetermined systems of PDEs. To find the complete 
symmetry group we will elaborate upon the Tanaka theory of symme- 
tries of vector distributions. 

We shall also discuss limitations of the theorem, by showing examples 
of degenerate systems when the external and internal symmetries are 
different. 



1 This is the parabolic subgroup P2 C G2 giving the contact gradation. 
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I. Symmetries of differential equations 

A differential equation £ is often treated geometrically as a subman- 
ifold in the space of jets J k (B, F), where B is the space of independent 
variables x = (x 1 , . . . ,x n ) and F is the space of dependent variables 
u = (it 1 , . . . , u m ), which we write as u = u(x). A choice of these coor- 
dinates on J° = B x F yields canonical coordinates on the space J k : 
(u J a ), where multiindex a — (zi, . . . , i n ) has length \a\ — J2a=i — ^- 

We will assume regularity, i.e. that the projection maps tt^i : J k — > 
J 1 have constant ranks when restricted to £ and its prolongations. 
Usually it is assumed that the equation has pure order k and then 
7Tfc,fc-i : £ — > J fc_1 is surjective. But we will allow PDEs in the system 
£ to have different orders (of which k is the maximal). For the setup 
of this theory we refer to |KL] . 

The Cartan distribution Ck on J k is given as the kernel of the forms 
9 J a = dv? a — Xir=i u i+iA xl f° r & U J e an d H < k. It is 

generated as follows: 

C k = (V x i = d x i + u i+iAi' d ui '■ 1 < ^ < W\ = k ) 

j;\r\<k 

Cartan distribution C\ on J 1 can be lifted to J k via n^]. In particular 
Ci is the classical contact distribution. 

These distributions for different / are related by the formula^ 

[Ci, C{\ = Ci-i. 

The classical Lie-Backlund theorem claims that a symmetry of dis- 
tribution Ck on J k is necessarily the lift of a contact transformation 
from J 1 if m = 1 or a point transformation of J° if m > 1. The proof 
follows almost immediately from the above display formula. 

A symmetry of (J k ,Ck), preserving £, is called an external symmetry 
of the system. The Lie algebra of all external symmetries is denoted 
by Sym(£). The internal symmetries Sym int (£) are by definition the 
symmetries of the induced distribution = Ck D T£ on £. 



For brevity sake we will write [V, A] instead of more appropriate "distribution 
whose module of sections equals [r(V),r(A)]" (regularity assumed). Notice that 
we always have [V, A] D V + A. 
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Restriction obviously gives the homomorphism Sym(£) — > Sym ini (£). 
Generalized Lie-Backlund theorems address the occasions when it is an 
isomorphism. 

In what follows we will consider an overdetermined system of PDEs 
£. It will be assumed formally integrable, meaning that for all prolon- 
gations £k+t = £® (given as the locus in J k+t of defining relations for £ 
differentiated < t times) the projections Ttk+t,k+t-i '■ £k+t — > £k+t-\ are 
afline bundles. In other words, all the compatibility conditions vanish. 

For technical reasons we will need to prolong to the level s = k + t, 
when £ s is involutive. This means vanishing of the Spencer cohomology 
H l, i{£) = 0, Vz > s,j > 0, or equivalently fulfilment of the Cartan test 
|C3[ IKLV] . Then for i > s the symbol 

Qi = Ker(c/7Tj i i_i : T£ { -> T£;_i) 

has dimension growth in accordance with Hilbert polynomial |KLj . If 

Xi G £i and n it0 (xi) = (x,y), we can identify gi{xj) C S l T*B ® T y F. 

2. Lie class u — 1 compatible systems 

The (complex) characteristic variety can be defined as projectiviza- 
tion of the set of complex characteristic covectors 

Char£(£) = F{p E T* X B ® C : p k ® T y F n ^ 0}. 

If Char c (£) = then the system is of finite type, i.e. its solution 
space if finite-dimensional. The next by complication case is when 
Char c (£) consists of one point counted with multiplicity (then the 
characteristic variety must be real). These are the systems of Lie class 
u = 1 and their solutions are parametrized 1 function of 1 argument!!. 

Another way to describe such systems is the following. For large 
% (precisely starting from the level when Si becomes involutive) the 
symbol of the system stabilizes: dim^j = 1. We refer to IK2I for a 
detailed discussion of these systems. 

By a theorem of S. Lie compatible PDE systems of class u = 1 
are integrable via ODEs [Lj . The proof given in |K 3 | uses the following 
observation: For such systems £, starting from the involutivity jet-level, 
the Cartan distribution Cs has rank equal to (n + 1) and it contains 
the subspace II of Cauchy characteristics of rank (n — 1). Since the 
prolongation, required to achieve involutivity, changes neither the class 
lo nor the external symmetry group, we will assume that already £ is 
involutive (this does not restrict generality of the results). 



In terminology of Elie Cartan uj is the Cartan integer s\ (provided the Cartan 
character is 1: S2 = 0). 
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The (local) quotient by the foliation tangent to II maps (£,Cs) to 
a manifold with rank 2 distribution. Let us call this pair (M, A) the 
reduction of our system. 

Shifts along Cauchy characteristic direction are trivial symmetries of 
£, so the space of internal symmetries is always infinite (while this for 
external is usually not). To compensate this we consider the symmetry 
algebra of the reduction Sym(M, A) (all our considerations are local). 
Since the subalgebra of trivial symmetries is an ideal in Sym int (£) , we 
have the homomorphism induced by the restriction 

Sym(£) — »■ Sym(M, A). (1) 

The goal of the paper is to demonstrate that this is an isomorphism 
(under certain "general position" assumptions). 

For simplicity of the exposition we restrict in the next section to the 
case of scalar PDEs (m = 1) on the plane (n = 2). We will denote a 
system of differential equations of orders fei, . . . , k r by Y7i 



A zoo of such systems is given in IK2I - In this reference and IK3I it 
is shown that internal geometry of linearizable systems is quite simple: 
their reduction correspond to Goursat distributions. 

Recall that for a distribution A the weak derived flag is defined by 
{Ai = A, A i+ i = [A, A,,]}, and the strong derived flag is given by the 
formula {Vi = A, Vj+j = [V;,V;]}. The sequences {dim(A i+1 /Aj)} 
and {dim(Vj + i/Vj)} will be called weak resp. strong growth vectors. 

Goursat distributions have both weak and strong growth vectors 
(2, 1, 1, . . . , 1) and are isomorphic to the Cartan distribution on the jet 
space J k (R, R). Thus the internal geometry of the linear class u = 1 
overdetermined compatible systems is trivial. 

On the other hand the external geometry (which is governed by 
the pseudogroup of point triangular transformations) is rich and is 
characterized by differential invariants. 

An interesting example of systems of type E 2 + E 3 can be found 
already in Cartan ( IC2I , p. 147). Modifying his PDEs a bit we get: 

£ ■ {^XIJ 0, Uyyy 0}, £ . \Uyy 0, U XX y } . 

Internally £ ~ £ since the growth vectors of both reductions are 
(2,1,1,1), but the systems are not equivalent externally because the 
second order equations are hyperbolic and parabolic respectively. Thus 
generalized Lie-Backlund theorem fails in this case. By |0 2 | such situ- 
ation is not possible for systems of type 2E 2 . 
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3. TWO EXAMPLES OF CLASS CO = I 

Let us study some partial cases with n = 2 independent variables. 
We will use the classical notations p = u x , q = u y , r = u xx , s = u xy , 

t ^yyi *y ^xxxi ft y>xxy-> T ^xyyi ^yyy' 

Consider at first the system E 2 + E 3 . Without loss of generality we 
can assume that the characteristic^ is d y — A d x , where A is a function 
on 2-jets (more precisely on equation E%). Thus the system £ writes 

t = F(r, s,...), (3 = A(r, s,...)a + B(r, s, . . . ), 

where dots mean terms of order 1 (and A satisfies A 2 = F S A + F r ). The 
Cauchy characteristic vector field for the distribution is £ = T> y — 
AT> X + gd a for some function g, and it is transversal to E : y = const. 

Consequently the reduced rank 2 distribution (on the quotient M) 
can be interpreted as the following distribution on E: A 2 = fl TE 
(the value of constant plays no role). In the canonical coordinates 

A = (V x = d x + pd u + r d p + s d q + ad r + (Aot + B) d s , d a ) . 

This distribution has growth vector (2, 3,4,...) and so is de-pro longable, 
i.e. locally A = P(A) for some rank 2 distribution A on 6-dimensional 
manifold M. Indeed, d a is the Cauchy characteristic for the derived 
distribution A 2 . 

Thus the weak derived flag {Aj} differs drastically from the strong 
derived flag {Vj}. One possibility is to perform de-prolongation (in 
this case it is easy - to restrict A 2 to the transversal M : a = const), 
but we will instead consider the strong derived flag. It has the following 
generators, provided the strong growth vector is (2, 1, 1, 2, 1)0 

Vi = (ei,d a ), ei = -V x 

V 2 = Vi + (e[), ei = [ex, d a ] = d r + Ad s 

V 3 = V 2 + (e 2 ), e 2 = [e 1; e[] = d p + Ad q + (e[(B) - V x (A))d s 

V 4 = V 3 + (e 3 , e'g), e 3 = [e u e 2 ] = d u + . . . , 

e' 3 = [e' 1 ,e 2 ] = e' 1 (A)d g + Cd s , 
V 5 = V 4 + (e 4 ) = TM, e 4 = [ei, e 3 ], [d, e' 3 ] = [e[, e 3 ] or [ei, e' 3 ], 

where the coefficients, like C or dots, will not be indicated explicitly 
(the growth vector of the reduction on M is (2, 1,2, 1) with the flag 
generated by vectors (ei, ei; e 2 ; e 3 , e' 3 ; e 4 )). 



4 This vector belongs to the annihilator of the characteristic covector dx + Ady. 
5 This is the case of general position; elsewise there are more de-prolongations or 
there exists a first integral of the distribution. 
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If we know the vertical distribution^] V 2 = (e[, d a ), we can construct 
the distribution V4 = [V 2 , V3] = (ex, e' x , e 2 , e' 3 , d a ). 

Distribution A will be called sufficiently non- degenerate if V 5 = 
[V' 2 ,V 4 ] has rank 6 (i.e. we can use [e^eg] for e 4 ). In this case it 
contains the fiber (d p , d q , d r , d s , d a ) of the projection G?7r 3j o : T£ — > TJ° 
as a codimension 1 subspace, so we re-cover the contact distribution 
on J 1 (lifted to £): & = V 5 + (£)■ 

The second particular type we want to examine is the system 3E 3 , 
which writes 

(3 = F(a,...), y = G(a,...), 5 = H(a,...), 

now dots mean terms of order 2. The condition of 1 common charac- 
teristic implies G a = F%, H a = and we assume £ fully nonlinear 
F<xa 7^ (this is sufficient non- degeneracy] the quasi-linear case yields 
formulae similar to the first example, so most of them also work). 

The characteristic is then d y — F a d x , whence the Cauchy character- 
istic of Cs is £ = V y — F a V x + q d a for some function g on £ and we 
again choose the transversal as S : y = const. 

This will be identified with the quotient M and the induced rank 2 
distribution is 

A = (V x = d x + pd u + r d p + s d q + ad r + Fd s + Gd t , d a ) . 

Now due to full nonlinearity the square of this distribution does not 
have Cauchy characteristics (i.e. it is not de-prolongable) . This implies 
that the first 3 elements of the weak and strong derived flags are the 
same, and so in our cas^ both growth vectors are (2, 1, 2, 3). 
The flags are given by 

Vi = (ei, e[), e x = -V x , e[ = d a 

V 2 = Vi + (e 2 >, e 2 = [ei, e[] = d r + F a d s + F 2 a d t 

V 3 = V 2 + (e 3 , e' 3 ), e 3 = [ei, e 2 ] = d p + F a d q + Cd s + Dd t , 

e' 3 = [e[, e 2 ]/F aa = d s + 2F a d t , 

V 4 = V 3 + (e 4 , e' 4) el) = TM, e 4 = [e x , e 3 ] = d u + . . . , 

< = [ ei , e' 3 ] = d q + C d s + D d t , e'i = [e[, e' 3 ]/(2F aa ) = d t . 

Thus if we know the vertical distribution (e^) C A we can re-cover the 
Cartan distribution on J 1 : C\ — [e'^Vs] + (£) (here we use the fact 
that e' A = [e[,e 3 ]/F aa modV 3 ). 

6 Notice that it is in the kernel of the bracket-map A 2 V2 — > V3 since d a is the 
Cauchy characteristic for V2. 

7 This equality of weak and strong derived flags can fail already when the length 
of the weak growth vector is 5, even in fully nonlinear case. 



8 



BORIS KRUGLIKOV 



Another way is to re-cover the contact distribution on J 2 : C 2 = 
K, [ei, V 2 ]] + (0 and then get d = [C 2 ,C 2 ]. 

Of course, in both cases we do not know the vertical distributions 
from internal viewpoint. However generically (if the system is suf- 
ficiently non-linear) they rotate inside the corresponding internally- 
canonical distributions when moving along the Cauchy characteristic. 
This crucial observation will allow us to prove the equivalence. 

4. Sufficient non-degeneracy 

We need to specify what conditions on £ should be imposed to obtain 
the generalized Lie-Backlund theorem. The first one is: 

(N) Nonlinearity of £ A is not Goursat. 

Thus the strong growth vector of A is (2, 1, . . . , 1, 2, . . . ), i.e. there 
exists s > 2 such that dim(V s /V s _i) = 2 and this is the first 2 after a 
sequence of 1 (s is minimal). 

Then dim(V s ) = s + 2 and this distribution has (s — 3)-dimensional 
sub-distribution C V s _3 of Cauchy characteristics for V s (that 

are also Cauchy characteristics for V s _ 2 ). 

Denote the projection along Cauchy characteristic space II C Cs by 
zu : £ — y M; it maps Cs — > A. The strong growth vector of Cs is 
(n + l,l,...,l,2,...), where the first dimension 2 occurs at the place 
s. The strong derived flag is V, = G7 _1 (Vj). Let O s -3 = ^ _1 (n s -3) C 
V s _3 be the space of Cauchy characteristics for V s (and V s _2). 

Let v s _2 = V s _2nKer(<i7r), where n : £ — y B is the natural projection 
in jets; rank(t; s _2) = s — 2. Then 

□ a -3 C v s -2 + n C v s _ 2 + A = V s _ 2 , (2) 

and both inclusions have codimension 1: rank(D s _3) = n + s — 4, 
rank(-L> s _ 2 + 11) = n + s — 3, rank(V s _2) = n + s — 2. 

Since II is the space of Cauchy characteristics, we have [IT, V s _ 2 ] = 
V s _2- Our second requirement is that the space -u s _ 2 rotates along II: 

(R) [n,t; s _ 2 ] = V s _ 2 . 

This condition holds if the system exhibits some degree of non-linearity; 
for instance, a system of type kEk with n = 2 satisfies (R) whenever 
the equations in it are not all quasi-linear. 

Denote the operation of taking bracket with v s -2 by ad„ s _ 2 . We 
have ad„ s _ 2 (V s _i) C V s . Let us continue taking brackets and denote 
the limit by ad!^_ 2 (V s _i). Our third requirement is that this latter 
generates the Cartan distribution: 

(G) (k — s+ l)-strong derived of ad^_ 2 (V s _i) is equal to C\. 
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Actually, if (k — s + 2) is less than the minimal order of £, we can 
formulate the condition in a simpler form: 

(G') ad£_ 2 (V s _i)=C fe _ s+2 . 

This implies, in particular, that the Cartan distribution Cs on £ is 
completely non-holonomic, i.e. its iterated brackets generate T£. 

Finally let us consider the condition on the space of Cauchy charac- 
teristics, which arise only in the case n > 2 (when dim II = n — 1 > 1). 
Since the inclusion II + v s ^2 C [II, V s _ 2 ] = V s _ 2 has codimension 
1, condition (R) implies existence of codimension 1 sub-distribution 
III C II such that [n 1 ,u s _ 2 ] C II + -u s _ 2 . We assume, in addition, to 
(jR), that we have equality instead of the general inclusion. 

Then there exists a codimension 1 sub-distribution II 2 C Hi such 
that [II 2 ,i; s _ 2 ] C IIi + v s _2- Again, we strengthen this to be equality 
and continue. To summarize we arrive to the filtration 

n = n d iii d n 2 • • • d n n _ 2 , 

where rank(ILj = n — i — 1 and [nj,u s _ 2 ] C IL;_i + v s -2- Our last 
condition, strengthening (R), is that the filtration rotates along v s -2, 
i.e. the last defining relation is the equality: 

(R+) [IT, v s - 2 ] = Ilj-i + V s -2, 

where we let II_i = Cg = ro _1 (A) to include the condition (R). Notice 
that there are no difference between (R) and (R+) for n = 2. 

Definition. An over determined compatible system £ of class u> = 1 is 
sufficiently non- degenerate if it satisfies conditions (N), (R+) and (G). 

Now we can prove the main Theorem. 

5. Proof of the main result 

Let us start by considering the case, when the base B has dimension 
n = 2, so that the space II of Cauchy characteristics is 1-dimensional. 

It is obvious that an external symmetry of £ induces an internal 
symmetry of the reduction (M, A). We have to show that the inverse 
exists and is unique. 

Since II is the space of Cauchy characteristics, we have (8,0s) ~ 
(M, A) x (II, II). This local diffeomorphism does not show the contact 
structure. The latter can be uncovered due to condition (R ) as follows. 

By (j5J) we can locally identify the space (8, v s -2 + II) with P s _ 2 = 
{V s ~ 2 : n s _3 C V s ~ 2 C V s _ 2 }. Since every (s — 2)-dimensional space 
ys-2 £ p s 2 j s squeezed between the subspaces of TM of dimensions 
(s — 3) and (s — 1), this P s - 2 is fibered over M with fibers of dimension 
1, and the condition (R) means that the line field V s ~ 2 /D s _3 rotates 
along the fiber, which can be locally identified with IT. 
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Indeed, the above construction corresponds to the prolongation of 
rank 2 distributions as follows. Recall that for a rank 2 distribution A 
on M its prolongation M is the manifold of all 1- dimensional subspaces 
£ C A, with the natural projection p : M — > M. The fiber of p over 
x G M is PA 5 ~ S* 1 and the natural lift of the distribution is given by 
the formula A e = dip~ l {f). This is a rank 2 distribution on M with 
the derived rank 3 distribution equal to [A, A] = dp~ l (K). The space 
of Cauchy characteristics of the latter is the fiber TS l . 

Now letting (M, A) = (M, V s _ 2 )/D s _ 3 = (£, V s _ 2 )/l%_ 3 (this is 
a local construction, quotient by D s _3 is possible as it is the space of 
Cauchy characteristics), we can locally identify the distribution -u s _ 2 +n 
on £ with the distribution A x D s _ 3 on M x R s ~ 3 . 

This construction allows us to uniquely (locally) recover the pair 
(V s -2,v s -2) on £ from the internal geometry of the distribution A on 
M. Since, by condition (G), the pair (V s _ 2 , f s - 2 ) generates the contact 
distribution on the space of jets, the symmetries of (M, A) are bijective 
with external Lie infinitesimal transformations of £. 

Consider now the case n > 2. By definition the bracket with IT acts 
trivially on the distribution V s _ 2 /O s -3. The previous construction 
identifies the line bundle IT/Tli with the projectivization of the rank 2 
distribution A = V s _ 2 /D s _3. Thus the distribution A + IT/IIi can be 
identified with the square [A, A] of the prolonged distribution A. This 
latter is identified with t; s _ 2 /D s _ 3 + IL/Hi, with IT/IIi corresponding 
to its vertical line subbundle. 

By condition (R+) this line bundle rotates along Ili/n 2 and does 
not rotate along Il 2 . Thus we can identify the line bundle Ili/n 2 with 
the prolongation of A. Continuing in the same way we identify IT with 
(n — l)-st iterated prolongation of the rank 2 distribution A over M, 
and so we uniquely recover the contact geometry of £ from the internal 
geometry of (M, A). QED. 



(we use jet multi-index notations u w = u x , u 01 = u y , w 20 = u xx etc; 
parameter A is to be excluded). 

For k = 1 this is the equation 2u y = u x , which is equivalent via a 
potential change of variables to the equation of gas dynamics u y = uu x . 
Both have infinite-dimensional contact symmetry algebra cont(lR 3 ). 



6. Application I 



Consider the following overdetermined system £ k C J fc (M 2 ): 
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For k = 2 this is the celebrated involutive second order PDE system 
considered by E.Cartain in 1893 and 1910 [Ui 



/ - - ^!\ 

He has shown that the contact symmetry group is the (split) excep- 
tional Lie group G 2 - 

Theorem 1. The algebra of contact symmetries of the equation S k for 
k > 2 has dimension k(k + l)/2 + Q and is isomorphic to the semi-direct 
product rife x gl 2; where n k is a nilpotent Lie algebra of length k. 

Every contact vector field on J 1 = J : (R 2 ) has the form (here T)\ = 
d x + u x d u , V\ = d y + Uyd u ): 

X f = -fuM ~ fuV\ + fd u + Vl{f) d Ux + via) d Uy , 

where / G C 00 ^ 1 ) is the generating function (contact Hamiltonian) . 
Let us indicate the generating functions / of a basis of Sym(£ fc ): 

x l y 3 (i+j< k), u x , u y , u + y- u y , (k + l)u -xu x , 

y ■ u x + x k /k\, (k — l)y ■ u — xy ■ u x — y 2 ■ u y — x k+1 / (k + 1)! . 

Since the prolongation of the field Xf to the space J k is given by the 
formula (a is the multi-index of the derivation, T> a is the iterated total 
derivative and T>\ = Yl u cT+ud Uc , is the truncated total derivative 

\a\<k 

[KLV]) 

x f = -UM - fu y v k y + v*(f)d u „, 

\a\<k 

we see that the first elements of the above collection act trivially on 
itk-i,i- The same concerns the translations u x ,u y . The next elements 
in the first line are scalings and it is trivial to check they preserve the 
system £ k . The last elements (second line) of the above collection have 
the following prolongations: 

— y 9 X + -jg- d u + (u_x)\ Ux ^ u v (fc-2)! ^ u xx u xx d Uxy + ^u xy d Uyy + 

+ d Uh + u kfi d Uk _ %1 + 2u k -i fl d Uk _ 2 2 + ■■■ + kux^-i d Uo k ; 
xyd x + y 2 d y + {{k - l)yu - ^f^y) d u + 

+ ((k - 2)yu x - |f) d Ux + ((k - l)u - xu x + (k - 3)yu y ) d Uy H 

- (x + yu k) o)d Uk o - ^ (i x Uk-i+i,i-i + (i + l)yu k -i,i) d Uk _. .. 

l<i<k 

Now it is straightforward to check these preserve the system E k . 
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Let us consider the reduction of S k by the Cauchy characteristic, 
which is the vector field V k — XV k + gd\, where A = u k> o and the 
precise value of g is not important. 

The quotient (M, A) is (locally) isomorphic to the intersection with 
the transversal {y = const}, so that A = (V k , d\), where 

V k = d x + u w d u + u 2 q d uw + Mn d U01 H h 

+ Xd Uk _ 10 + \\ 2 d Uk _ 21 H h |A d^^. 

Denoting w\ = Uij, < % + j < A; — 1, we see that the above rank 2 
distribution corresponds to the following underdetermined ODE system 
on w = (w°, . . . , as a function of x (we understand = d l x {w^))\ 

y k : K = A, wU = ±A 2 , . . . , wt 1 = l k x k }. 

By our main theorem, the algebra of external symmetries Sjm cont (Sk) 
is isomorphic to the algebra of internal symmetries of the above Monge 
equation 34 C J k,k ~ 1 ''''' 2,1 (R, IR fc ) (considered as a submanifold in the 
mixed-jet space jl, which is identical with the algebra Sym(A). 

The latter has the following basis (we write only the point part of 
the transformation, from which it can be uniquely recovered by the 
prolongation): 

x = d x , wi = \x l d wJ (o<i + j<k), 

L = ^x k d w o + d w i + 2w\ d w 2 H h (k - l)w k ~ 2 d w k-i, 

R = j^x k+1 d w o + (xw? - (k - l)w°) d w i + 2(xwl - {k - 2)w l ) d w 2+ 

+3(xw 2 -(k- 3)w 2 ) d w 3 H h (k - l)(xw k ' 2 - w k ~ 2 ) d w k-i, 

S x = x d x + kw° d w o + (k - l)^ 1 d w i H h 2w k ~ 2 d w k-2 + w fc_1 d w k-i, 

S 2 = w° d w a + 2w x d w i + 3w 2 d w 2 + h kw k ~ x d w k-i , 

T = A d x + (Aw? - w l ) d w o + (\w{ - w 2 ) d w i+ 

+ ■■■ + {Xwt 2 - w k ~ l ) d wk -2 + fe^A fc+1 d wk -i . 

Notice that (X, W- , L) is a translational part and (R, Si, S 2 , T) is the 
stabilizer of the origo (all coordinates vanish) o G M (to check this 
for L one has to prolong the formulae to see the term +d\; recall that 



8 Notice that y 2 C J 2a (R,R x E) is the Hilbert-Cartan equation. In [XK] it was 
generalized in some aspects, different from the present paper. 
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A = Wfr). The non-trivial commutators are the following 

[x, l\ = w£_ lt [x, w}\ = wu [l, w(\ = -(j + (i > 0), 

[X, R] = L, [Wj, R] = (j + l)(i + j + 1 - 

\wi s 1 ] = (k-i- mi wi sy = (j + & sy = a 

[L,T] = X, [W^T] = -wr l (j > 0), [R,T] = St - 5 2 , 
[R, Si] = -R, [R, S 2 ] = R, [T, S^ = T, [T, S 2 ] = -T. 

This algebra is graded with the following weights: w(X) = w(L) — — 1, 
w{Wi) —i — k — 1, w{R) = w{T) = w(Si) = w{S 2 ) = 0. 

Even more, it is bi-graded: b{X) = (-1,0), b(L) = (0, -1), 6(Wf) = 
{i + j-k,-j-l), b(R) = (1,-1), 6(T) = (-1,1), 6(5-0 = (0,0), 
6(5*2) = (0, 0); the grading w is the total weight of b. 

It is rather straightforward to check that the above fields are sym- 
metries of the distribution A. This gives the lower bound for Sym(A). 
To get the upper bound one should calculate the Tanaka algebra of A. 

Recall [Tj that with every distribution one associates the sheaf of 
graded nilpotent Lie algebras (called symbol or Carnot algebras) m = 
®i<o9i, 9i = A_j/A_j_!, where {A^} is the weak derived flag of A. 
The point-wise bracket of m is induced by the commutator of vector 
fields, so that for every point x € M we get the Lie algebra m = m^. 

The Tanaka prolongation m = ©jj, is defined as the maximal graded 
Lie algebra with the given negative part m [Tj. It can also be defined 
via graded Lie algebra cohomology as Qi = Hl(m,m © Q ■ ■ ■ © 0i_i) 
for i > 0. In other words, Qi are constructed successively as maximal 
subspaces such that all possible Jacobi identities hold, see |Y 2 | . 

To calculate m one takes the generators T) k x , d\ of A, and computes all 
possible brackets. The resulting Carnot algebra is a graded nilpotent 
Lie algebra isomorphic to (X, L, W-). 

It can be described as a truncated double-graded free Lie algebra 
with fundamental part of grading —1 and rank 2. In the appendix we 
demonstrate that the Tanaka prolongation of is trivial in the positive 
grading. Hence the Tanaka algebra is n& © jjo, where q = gl(g_i). 

By the results of [H |Ki| and the Theorem from Introduction this 



gives the upper bound for the symmetries of both equations: 8k (PDE 
system) and 34 (ODE system). Since it coincides with the algebra 
of symmetries we already constructed, our description of Sym(^) ~ 
Sym(A) is complete. Theorem [1] is proved. 

Remark. We have demonstrated that A is the most symmetric distri- 
bution with the symbol m equal to the truncated double-graded free Lie 
algebra n^. By Theorem 4 of [AKj this implies that A is Tanaka-flat. 
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We can describe very symmetric systems, for which the correspond- 
ing distributions are non-flat. Let us start with the PDE models 
T h c J k (R 2 ): 

7 k = \ u k -i i = — : < i < k 

L im + 1 

which coincides with £ k for m — 1. For k = 2 this is (equivalent to) 
the family of involutive 2nd order PDE systems 

\ Hxt, A, U X y ■ ~, ^yy 



m + V yy 2m + 1 

with 7-dimensional contact symmetry algebra described by Cartan in 
|C 2 | (the next large after 14-dimensional G 2 ), see also |K 4 | . 
By the calculation similar to the above we get the following 

Theorem 2. The algebra of contact symmetries of the equation T k for 
k > 2 and generic m has dimension k(k + l)/2 + 4. Its basis consists 
of contact vector fields Xf with generating functions f as follows 

x l yi (0 < % + j < k), u x , u y , (km + l)u— mxu x , u + my u y . 

By the main theorem we can also represent this algebra as the inter- 
nal symmetry of the following underdetermined ODE system obtained 
from via reduction by the Cauchy characteristic: 

< = ^ wU = ^a-+\ . . . , «,*-* = v-^r^ k - 1)m+1 - 

(here = w^{x) are the dependent variables and w\ = d l x {w^)). 

The symmetry algebra has a filtration with the corresponding graded 
Lie algebra being the semi-direct product rife x M 2 , where rife is the same 
as in Theorem [TJ This associated graded algebra is a subalgebra of the 
Tanaka algebra fife of the flat model, IR 2 being the diagonal part of go- 



It is possible to show by the methods of K4 that the corresponding 



distribution A is sub-maximal symmetric with the Tanaka algebra fife. 



7. Application II 

The PDE system E k has solution space Sol(£fe) that is parametrized 
by 1 function of 1 argument and dim M — 2 = fc ^ ^ constants (these 
are the so-called Lie class u = 1 systems [Lj |K 2 [ IK3I ) . 

More general systems can be treated with the proposed technique 
too. Consider, for instance, the following system VJ^ C J fe (M 2 ), m < k, 
for which the right-hand side is the m-th tangent cone (of dim = m + 1 
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in of the normal projective curve defining Sk = 1Z° k : 

min(m,i+l) , ( 

K? = + £ ^-^A^ :0 <i<*} 

(the PDE system is obtained by excluding the additional parameters 
A, d, . . . , ( m and obtaining relations on the jet- variables Uk-i,i)- 

For instance for k = 3, m = 2 the system 1Z\ looks so (this is one 
highly nonlinear PDE of the 3rd order) 

A 2 

«« BB = A + Ci, «»v = T +ACi + C a , 

A 3 A 4 

U X yy = — + A'Cl + 2AC 2 , Uyyy = ~ + A^l + 3X^2. 

The system 1Z\ is the famous Goursat parabolic PDE on the plane; its 
contact symmetry group is G 2 the same as for (J3j). 

The system 1Z™ is involutive and its characteristic variety (both com- 
plex and real) consists of 1 point with multiplicity m + 1. The solution 
space Sol(7?.™) is parametrized by u> = m + 1 functions of 1 argument 
(and some constants). There is a reduction of this parabolic system 
via a characteristic involutive distribution (no longer a space of Cauchy 
characteristics) to a distribution A on the quotient that de-prolongs to 
a rank 2 distribution A (Monge system); this generalizes [G] . 

This reduction makes a bijection between contact symmetries of the 
PDE system and the ordinary symmetries of A. While this idea is 
more general, we will indicate only how it applies to our system 

Theorem 3. The algebra of contact symmetries of the equation 1Z™ 
for k > 2, < m < k, has dimension k{k + l)/2 + 6 and is isomorphic 
to the same Lie algebra x gl 2 as in Theorem [B 

Proof. Let us first show how to associate a rank 2 distribution to such 
a parabolic system of PDEs. 

Consider the Cartan distribution C of 1Z™. It is generated by two 
truncated total derivatives Vy and the vertical fields d\, d^, . . . , d^ m . 
The sub-distribution U m = (d^, . . . ,d^ m ) is integrable and we would 
like to quotient by it. 

It however does not commute with C, so we need to add the commu- 
tators 

Bd£(C)=C+(TH,...,rhn), 

where 771 = ^to ^ 9 u k -i ,v Vi+p = d\(Vi)- 

The quotient manifold is M = R fc ( fc+1 )/ 2+3 (x, y, A, u a : \a\ < k) and 
it is equipped with the distribution A = ad^ (C)/H generated by the 



'This R fc+1 (u fc o,...,w fe) = S k (R 2 )* is the fiber of 7r fe , fe _i : J fc (R 2 ) -> J^^R 2 ). 
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truncated total derivatives evaluated at Ci = • • • = Cm = 0, denoted 
V%., T3yi and by the fields rji, . . . , r) m , d\. For instance, for k = 3, m = 2 

V x = d x +pd u + rd p + sd q + Xd r + \\ 2 d s + §A 3 <9 t , 
V y = d y + q d u + s d p + 1 d q + \X 2 d r + §A 3 <9 S + \\%, 
rn = d r + \d s + \ 2 d t , r) 2 = d s + 2\d t 

The vector £ = T) y — XD X is the Cauchy characteristic of the distri- 
bution A. Let (M, A + ) = (M, A)/£ be the local quotient. It is easy to 
see that A + is the ra-th derived (both weak and strong) distribution 
of the rank 2 distribution A = (T> x , d\). 

Thus any contact symmetry of C descends to an ordinary symmetry 
of A. Conversely, Sym(A) = Sym(A + ) because prolongation preserves 
the symmetries, and then by our main Theorem the latter coincide 
with the contact symmetries of A. This distribution corresponds to 
the equation = TZ^, from which by taking the m-th tangent cone we 
obtain our system W£. Thus any ordinary symmetry of A uniquely 
lifts to a contact symmetry of C. 

The claim now follows from Theorem [T] on the symmetries of □ 

8. The case of n > 2 independent variables 

It is easy to see that the generalized Lie-Backlund theorem, as it 
was stated in Introduction, fails if we allow equations of the first order 
without non-degeneracy assumptions. Indeed, this is so with any class 
u — 1 system of the type IE2 + E\ in n > 2 independent variables. 

For instance, we can take ([3]) as 2E 2 and trivial E\\ 

P ^3 

U xx A, U X y ^ j V"yy g > ^2 

Then z H- Z(x,y, z,u,u x ,u y ,u z ) is a contact symmetry, inducing the 
trivial transformation of the reduction (M 5 , Ahc)> which is the Hilbert- 
Cartan equation. In fact, the general contact symmetry Xf has gen- 
erating function / = / + u z ■ f ', where fo = fo(x,y,u,u x ,u y ) is the 
14-parametric generating function corresponding to GVaction on ([3]), 
and / is an arbitrary function on J 1 (IR 3 ). Thus the generalized Lie- 
Backlund theorem fails: the map is not injective. 
Let's consider a similar system of the second order: 

A 3 

U'XX A, U x y ^ ) g ) 0, Uy z 0, U Z Z 0' 

Its reduction (along rank 2 distribution II) is the following rank 2 dis- 
tribution in 6D: (M,_A) = (M 5 , A HC ) x (R, 0). Let t be the first integral 
of the distribution A (it is not completely non-holonomic). The sym- 
metries are t 1— y T(t) and (^-transformations of the first factor, with 



LIE-BACKLUND THEOREM FOR LIE CLASS ui = 1 SYSTEMS 



17 



coefficients parametrized as functions of t. Thus in total Sym(M, A) 
is parametrized by 15 functions of 1 variable. 

A direct calculation shows that the contact algebra is precisely the 
same: (Lie(£?2) © M) M , so the generalized Lie-Backlund theorem holds. 

Similarly this theorem holds for higher order equations. For instance, 
the system of type 3E 3 + 3E 2 



A 2 A3 

U XXX - A, U XX y - 2 , U X yy - g , Uyyy - ^,11^-0, Uy Z - 0, U ZZ - 

reduces to (M 9 , A 2 ) = (Mi, Aj) x (R, 0), and both contact symmetry 
of the PDE system and the internal symmetry of the reduction are 
parametrized by 13 functions of 1 variable. More precisely, these alge- 
bras are both equal to (g©R) R , where g is the algebra of symmetries of 
(Mi, Ai) equivalent to the Monge system of 2 equations on 3 unknowns 
(the indices denote the number of x-derivatives) studied in Section [6j 

v 3 = A, w 2 = |A 2 , z\ = f A 3 . 

On the other hand if we consider the system of the 3rd order 9E 3 

A 2 A 3 _ A4 

Uxxx Uxxy 2 ' ^ x yy g ' ^yyy ^ > 

l^xxz 0; l^xyz 0; ^yyz ^' ^xzz 0, Uy Z z 0, U zzz 0, 

its contact symmetry algebra is 21-dimensional, while the internal sym- 
metry algebra of the reduction (the distribution there is not completely 
holonomic) is infinite-dimensional. Thus the generalized Lie-Backlund 
theorem fails: the map ([1]) is not surjective. 

All these systems have a kind of degeneracy, so let us consider a 
totally non-linear (and non-degenerate) system, for which our Theorem 
guarantees that the generalized Lie-Backlund theorem holds. 

We start with the equations of pure order 2 and Lie class u> = 1, for 
which the result, that contact external symmetries of the PDE system 
and the internal symmetries of the reduction coincide, follows also from 



Yi| . Consider the following example. 

C \mi+mj+l 

£ = \ ua = : 1 < i < j < n>. 

\ 3 rrii + mj + 1 - -J - j 

This system is involutive and, if m, ^ rrij for i ^ j, it has no first 
integrals. We can always achieve mi = by re-parametrization. The 
reduction of £ is given by 

U«2 + l ( n .X \m a +l („,1 \m n +l 

2 _ \ u xx) 3 _ \ u xx) n _ \ u xx) / a\ 

m 2 + 1 m 3 + 1 m n + 1 

(the superscript numbers the unknown functions). The internal sym- 
metry algebra is maximal for m 2 = 1, m 3 = 2, . . . , m n = n — 1, and 
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has dimension (2n + 5); the proof of this fact is similar to the proof 
of Theorem [T] and so is omitted. In fact, the rank 2 distribution of 
ODE system (J4]) with the prescribed parameters is internally equiva- 
lent to the most symmetric Monge equation [AK] (subscripts denote 
the derivatives) 



(its symmetry algebra with notation ti jTl = p„ + 3 was studied in [AKj ). 

The symmetry algebra is the same for the ODE and PDE systems for 
all rrii, but for generic parameters the algebra has smaller dimension. 

Now the above PDE system 8 can be generalized to higher orders. 
For the third order it writes as 



It is involutive, but to achieve complete non-holonomy (no first inte- 
grals), we have to assume that all the numbers rrij+mk, 1 < j < k < n, 
are different. The reduction can be again easily described. 

Consider a particular case n = 3, mi = 0, m<i = 1, = 4, satisfying 
the above restriction. The symmetry algebra of both ODE reduction 
and the PDE system is 15-dimensional. 

Similarly for n = 4 independent variables, and parameters mi = 0, 
m-2 = 1, m 3 = 4, m 4 = 15, we calculate both symmetry algebras to 
have dimension 21. 

We conclude that the generalized Lie-Backlund theorem holds in 
many interesting cases. The calculations are easier for the reduced 
ODE models, where Tanaka theory helps to calculate symmetries by 
using the algebraic methods. 

Appendix A. Truncated double-graded free Lie algebra 

A double-graded free Lie algebra is such a graded Lie algebra rioo = 
©i<o0« that its fundamental space is a direct sum of two subspaces 
IIi, n 2 , and its commutators freely generate the whole algebra with no 
other relations than the commutation of ad^ , ad„ 2 for v t G 11, (whence 
Z-grading can be refined to Z © Z-grading) . 

Let us specify this only in the case of current interest dimg_! = 2, 
when the fundamental space has basis eio, eoi according to the above 
splitting. Then g_ 2 is generated by e u = [ei ,e i], g_ 3 by e 2 i = 
[eio, en], ei2 = [eoi, en] etc. The only relations this infinite-dimensional 
Lie algebra rioo admits are [eio,e ij+ i] = [e i, ei+ij) = ej+ij+i. 

Thus 0_fc_i = (ejfe,i, efc-i^, . . . , ei t fc) has dimension k. It's easy to 
check that rioo is a Lie algebra. Many interesting graded nilpotent Lie 
algebras are quotients of this algebra (see some in [AKj ). 



V\ = {Znf 
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The truncated algebra is = g_fc © • • • © (the brackets between 
0_i and g_j are zero if i + j > k). Our goal is to calculate the Tanaka 
prolongation of (beware: this sub-script k is not a grading). 

Theorem 4. For k > 3 the Tanaka prolongation of n k is supported in 
non-positive grading: = © q , g = gl(g_i). 

We consider prolongation in the usual sense, not preserving double- 
grading (in particular go consists of grading zero derivations, not nec- 
essary preserving the splitting of 0-i). 

Proof. We first state that g = gl(g_i) is the maximal possible algebra 
of grading preserving derivations. The easiest way to see this is to 
calculate the prolongation of an element 

a b 
c d 



h 



e 0o 



(the matrix in the basis eio,eoi) and check that as derivation it pre- 
serves the relations [ad ei0 , ad eoi ] = 0. 

This readily follows from the following formula that can be proved 
by induction (1 < i < k — 1): 

h{e k - i: i) = (i-l)b e fc _ i+M _i + ((£;-i)a-M d) e k - iti +(k-i-l)ce k -i-i d+1 . 

Let now u G gi. Since is fundamental (generates n^), this ele- 
ment is uniquely determined by specifying uj(e w ) = h', cu(e i) = h", 
where we denote = [ou,£\ and the elements of g have the form 



ti 



1/ 
d 1 



h" 



b" 
d" 



By the Leibniz rule we calculate: 

w(en) = {b' - a") e 10 + (df - c") e 01 , 

u(e 21 ) = (a' + 2d' - c") e n , u(e 12 ) = (2a" + d" - b') e n , 

w(e 3 i) = (3a' + 3d' - c") e 2 i + c'e 12 , 

w(e 22 ) = (2a" + d") e 2 i + (a + 2d') e 12 , 

w(eis) = b"e 21 + (3a" + 3d" - b') e 12 . 

In calculation of u(e 22 ) we can use two representations [ei ,ei 2 ] = 
[eoi, e 2 i] and this gives the same result. Thus if we truncate on the level 
k = 3 (so that e^i = e 22 = ei3 = 0), then all coefficients to the right in 
the last three lines vanish and we obtain that q± is 2-dimensional (in 
agreement with the known grading of the exceptional Lie group G 2 ). 
Let now k > 3 and we calculate the action on $j_ 5 . We get 

w(e 4 i) = uj([e 10 , e 31 ]) = (6a' + Ad' - c") e 31 + 3c'e 22 , 

u(e u ) = w([e i, ei 3 ]) = 36"e 22 + (4a" + 6d" - b') e 13 
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uniquely and 

w(e 32 ) = w([ei , e 22 ]) = (2a" + d" + 6') e 31 + (3a' + 4d') e 22 + c'e 13 
= w([e i, e 3 i]) = (3a" + d") e 31 + (3a' + 3d' + c") e 22 + c'ei 3 , 

w(e 23 ) = w([ei , e 13 ]) = 6"e 31 + (3a" + 3d" + 6') e 22 + (a' + 3d') e 13 

= w([e i, e 22 ]) = 6"e 31 + (4a" + 3d") e 22 + (a' + 2d' + c") e 13 

non- uniquely, which implies b' = a", c" = a". 

Further calculations show that the derivation respects the higher 
commutation relations, and we obtain by induction: 

, . iU — l).., /,, ... .. + ,,A 
w(e fe _ M+ i) = — b efc-i+i^-i + [{k - %)% a H d J e fc _ M 

+ ( (fe ~ i)(fe 2 ~'~ 1} a' + (< + - < - l)d') 

(fc-.-l)(A;-i-2) 
H c efc_j_ 2) j +2 . 

Now it's obvious that truncation on the level k, i.e. letting ek+i-i,i = 
for 1 < % < k, forces hi = h" = 0, so that g 1 = 0. □ 
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